r 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  074-0188 

Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other' aspect  of  this  collection  of  information,  including 
suggestions  for  reducing  this  burden  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302, 
and  to  the  Office  of  Management  and  Budget.  Paperwork  Reduction  Project  (0704-0188).  Washington.  DC  20503 

1 .  AGENCY  USE  ONLY  (Leave  blank)  2.  REPORT 

DATE  1997 

3.  REPORT  TYPE  AND  DATES  COVERED 

Scientific  paper 

4.  TITLE  AND  SUBTITLE 

Application  of  the  GRI  1.2  Methane  Oxidation  Model  to  Methane  and  Methanol  Oxidation  in 
Supercritical  Water 

5.  FUNDING  NUMBERS 

N/A 

6.  AUTHOR(S) 

Steven  F.  Rice 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Combustion  Research  Facility 

Sandia  National  Laborartories,  MS-9052 

P.O.  Box  969 

Livermore,  CA  94551-0969 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

N/A 

9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESSEES) 

SERDP 

901  North  Stuart  St.  Suite  303 

Arlington,  VA  22203 

10.  SPONSORING  /  MONITORING 
AGENCY  REPORT  NUMBER 

N/A 

11.  SUPPLEMENTARY  NOTES 

Proceedings  of  the  Fourth  International  Symposium  on  Supercritical  Fluids,  May  1997.  This  work  was  supported  in  part  by  SERDP. 
The  United  States  Government  has  a  royalty-free  license  throughout  the  world  in  all  copyrightable  material  contained  herein.  All 
other  rights  are  reserved  by  the  copyright  owner. 

12a.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Approved  for  public  release:  distribution  is  unlimited  | 

380710 

022 

1  12b.  DISTRIBUTION 

CODE 

A 

13.  ABSTRACT  (Maximum  200  Words) 

The  Gas  Research  Institute  (GRI)  has  been  leading  an  effort  over  the  past  few  years  to  consolidate  recent  developments  in  the 
elementary  reaction  modeling  of  the  oxidation  of  methane  for  combustion  applications  into  a  single  optimized  set  of  reactions  within 
a  standard  computational  framework.  Our  interest  in  the  translation  of  combustion-based  models  to  oxidation  in  supercritical  water 
has  led  us  to  examine  the  applicability  of  this  mechanism  to  the  low-temperature  high-density  conditions  of  supercritical  water 
oxidation  at  25  MPa  and  400-650  °C.  This  paper  shows  that  the  GRI  1.2  mechanism  accurately  represents  the  available  experimental 
results  on  methane  over  a  wide  temperature  and  concentration  range.  The  oxidation  of  methanol  is  not  well  represented  by  the  GRI 
mechanism  when  left  unchanged.  However,  if  a  key  modification  is  made  to  the  reactivity  of  H02,  good  agreement  with  the  methanol 
oxidation  results  is  achieved.  Although  designed  for  conventional  combustion  conditions,  the  GRI  mechanism  can  be  successfully 
extended  with  very  little  modification  to  much  lower  temperature  and  extreme  pressure  conditions. 

14.  SUBJECT  TERMS 

Supercritical  Water,  Oxidation,  Methane,  Combustion,  SERDP 

15.  NUMBER  OF  PAGES 

4 

16.  PRICE  CODE  N/A 

17.  SECURITY  CLASSIFICATION 

OF  REPORT 

unclass. 

18.  SECURITY 
CLASSIFICATION 
OF  THIS  PAGE 

unclass. 

19.  SECURITY  CLASSIFICATION 

OF  ABSTRACT 

unclass. 

20.  LIMITATION  OF 
ABSTRACT 

UL 

NSN  7540-01-280-5500  Standard  Form  298  (Rev.  2-89) 


Prescribed  by  ANSI  Std.  Z39-18 
298-102 


Lpnc  quality  inspected  i 


4th  International  Symposium  on  Supercritical  Fluids 

Application  of  the  GRI  1.2  Methane  Oxidation  Model  to  Methane  and 
Methanol  Oxidation  in  Supercritical  Water 

16-1997 

Steven  F.  Rice 

Combustion  Research  Facility 
Sandia  National  Laboratories,  MS-9052 
P.O.  Box  969 

Livermore,  CA  94551-0969 
Fax:  (510)  294-1004  Email:  sfrice@sandia.gov 

The  Gas  Research  Institute  (GRI)  has  been  leading  an  effort  over  the  past  few  years  to  consolidate  recent 
developments  in  the  elementary  reaction  modeling  of  the  oxidation  of  methane  for  combustion  applications  into 
a  single  optimized  set  of  reactions  within  a  standard  computational  framework.  Our  interest  in  the  translation  of 
combustion-based  models  to  oxidation  in  supercritical  water  has  led  us  to  examine  the  applicability  of  this 
mechanism  to  the  low-temperature  high-density  conditions  of  supercritical  water  oxidation  at  25  MPa  and  400  - 
650  °C.  This  paper  shows  that  the  GRI  1.2  mechanism  accurately  represents  the  available  experimental  results 
on  methane  over  a  wide  temperature  and  concentration  range.  The  oxidation  of  methanol  is  not  well  represented 
by  the  GRI  mechanism  when  left  unchanged.  However,  if  a  key  modification  is  made  to  the  reactivity  of  H02, 
good  agreement  with  the  methanol  oxidation  results  is  achieved.  Although  designed  for  conventional 
combustion  conditions,  the  GRI  mechanism  can  be  successfully  extended  with  very  little  modification  to  much 
lower  temperature  and  extreme  pressure  conditions. 


Introduction 

Recently,  there  have  been  key  developments  in 
reactor  design  using  transpiration-wall  technology 
that  mitigate  corrosion  and  scaling  issues  for 
applications  of  supercritical  water  oxidation 
(SCWO)  to  a  variety  of  problematic  wastes  [1,2]. 
However,  it  is  likely  that  a  thorough  understanding 
of  the  operative  chemical  kinetics  for  reaction 
initiation  will  be  required  to  underpin  future 
applications  of  this  novel  design  strategy.  Continued 
improvement  in  reactor  design  for  energy-  and  cost- 
efficient  operation,  and  to  assure  safety  will  be 
hampered  without  better  predictive  models  for  the 
time,  temperature,  pressure,  and  concentration 
dependence  of  the  oxidation  process. 

Much  of  the  current  understanding  of  the  rates 
and  mechanisms  of  reactions  in  supercritical  water  is 
limited  to  a  small,  but  growing,  list  of  empirical 
global  expressions  for  simple  chemicals.  However, 
global  expressions  are  of  limited  use  in  the 
formulation  of  a  predictive  model.  To  be  valuable  as 
a  predictive  design  tool,  a  model  must  be  based  on  at 
least  a  quantitative  mechanism  incorporating  the  key 
elementary  reaction  steps. 

The  first  attempt  to  use  an  elementary  reaction 
model  to  describe  the  kinetics  of  oxidation  of  the 
simple  organics  in  supercritical  water  was  reported 
by  Webley  and  Tester  [3].  The  approach  was  to 
modify  a  successful  mechanistic  scheme  generated 
for  high-temperature  gas-phase  oxidation  for  high- 
density  effects  on  unimolecular  reaction  rates.  This 
mechanism  successfully  reproduced  the  existing 
methanol  oxidation  data,  but  failed  to  even 
approximate  the  experimental  results  that  had  been 


obtained  for  methane.  Recently,  others  have 
developed  more  extensive  elementary  reaction 
models  for  the  oxidation  of  methane  and  methanol  in 
supercritical  water  [4-7].  Similar  to  the  approach  in 
Ref.  3,  all  of  these  newer  mechanisms  have 
incorporated  direct  customization  of  several  key 
reactions  for  high  pressure. 

In  this  paper,  we  examine  the  application  of  the 
recently  developed  GRI  1.2  mechanism  [8]  to  the 
oxidation  of  methane  and  methanol  at  approximately 
25  MPa  and  420  -  630  °C.  This  mechanism  has 
high-pressure  considerations  explicitly  included. 
The  goal  of  this  paper  is  to  make  a  connection 
between  the  high-pressure,  low-temperature  results 
from  recent  experiments  to  the  current  best  effort  at 
representing  oxidation  of  Cl  species  at  combustion 
conditions.  We  have  adopted  an  approach  that 
establishes  this  relationship  with  the  smallest  amount 
of  customization  as  possible,  with  all  of  the 
parameters  for  both  the  mechanism  and 
thermodynamics  originating  from  a  single  well- 
documented  source,  and  the  computation  conducted 
within  the  framework  of  a  well-established  code. 

Methods 

The  calculations  presented  here  are  conducted 
within  the  Chemkin  II  computational  framework  [9]. 
The  GRI  1.2  mechanism  is  designed  for  Chemkin  II 
and  is  accompanied  by  a  Chemkin-compatible 
thermodynamic  data  base.  The  several  calculations 
reported  here  using  the  mechanism  of  Schmitt  [5] 
use  the  GRI  thermodynamic  data  base  with  the 
addition  of  properties  for  CH3O2,  and  several  other 
less  important  alkyl  peroxides,  calculated  using 
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THERM  [10].  These  species  are  included  in  the 
Schmitt  mechanism,  but  not  in  the  GRI  data  base. 

Oxygen  is  used  as  the  oxidizer  at  an  equivalence 
ratio  of  0.75  for  calculations  where  the  stoichiometry 
is  not  specifically  listed.  Moderately  lean  conditions 
such  as  these  are  likely  to  be  used  in  industrial 
SCWO  systems. 

In  some  instances,  comparison  with  experiment 
is  done  using  an  effective  first-order  rate  constant, 
keff,  defined  as  the  reciprocal  of  the  time  at  which 
the  initial  fuel  mole  fraction  has  fallen  to  1/e  of  is 
original  amount  from  the  expression 

-d[CH4]/dt=keff[CH4]  (1) 

Thus,  much  of  the  emphasis  in  this  paper  is  on 
oxidation  rates  at  relatively  low  conversion. 
Unfortunately,  this  paper  simultaneously  illustrates 
that  these  reactions  generally  are  not  first  order  when 
the  fuel  concentration  is  varied  over  a  wide  range. 
Despite  this  problem,  keff  is  still  the  simplest  way  to 
represent  reaction  rates  when  the  reaction  conditions 
being  considered  result  in  the  rates  varying  by  many 
orders  of  magnitude. 

Results  and  discussion 

Fig.l  shows  a  comparison  of  the  results  from 
calculations  of  the  oxidation  rate  of  methane  with 
the  presently  available  experimental  data.  The 
results  for  390  -  440  °C  are  calculated  from  a  global 
fit  to  dozens  of  experiments  that  were  conducted  at 
27  MPa  and  an  initial  fuel  concentration  of 
0.1  mol/1  [1 1],  A  global  fit  to  these  data  produces  the 
relationship 

-d[CH4]/dt=  k  [CH4]l-84[02]-006  (2) 


and  k  =  1011  ‘expC-llSOO/T).  However,  Webley  and 
Tester  also  provide  a  calculation  for  a  pseudo  first- 
order  rate  constant,  calculated  form  the  same  data. 
The  fitted  values  for  that  expression  are  used  here 

for  W 

Temperature  (°C) 


Figure  1.  Comparison  of  calculated  effective  first 
order  rate  constants  for  the  oxidation  of 
methane  from  the  GRI  1.2  mechanism 
and  experimental  results. 

The  agreement  of  the  Chemkin  calculation  using 
the  unmodified  GRI  mechanism  with  the 
experimental  data  is  excellent.  Although  the 
calculation  reveals  fuel  consumption  curves  that 
exhibit  a  brief  induction  period,  it  is  a  small  fraction 
of  the  1/e  time  constant.  The  GRI  1.2  mechanism 
can  accurately  reproduce  the  available  fuel 
consumption  rates  for  the  oxidation  of  methane  over 
the  entire  temperature  range  likely  to  be  used  in 
SCWO  applications. 


where  the  concentrations  are  in  mol/1, 
k  =  lO^^expC-SOlOO/T),  and  the  units  of  the 
preexponential  factor  are  chosen  so  the  rate  is  in 
mol/l-s.  The  pseudo  first-order  rate  constant,  keff,  is 
then  calculated  from  these  parameters  assuming  the 
dependence  on  CH4  concentration  is  unity  and  the 
dependence  on  oxygen  is  zero  at  an  initial  fuel 
concentration  of  0.003  mol/1.  This  procedure  yields 
a  keff  for  a  comparable  concentration  to  that  used  in 
Webley  and  Tester  [3].  The  low-temperature 
experiments  covered  both  lean  and  rich  conditions, 
but  no  oxygen  concentration  dependence  was 
observed.  In  addition,  no  significant  induction 
period  was  observed.  A  global  fit  to  the  high 
temperature  data  [3]  gives  an  expression  for  the  fuel 
consumption  to  be 

-d[CH4]/dt=  k  [CH4]°"[02]0-66  (3) 
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Figure  2.  Flux  diagram  of  the  major  carbon  species 
in  the  oxidation  of  CH4  at  500  °C  and 
25.0  MPa  in  supercritical  water  predicted 
by  the  GRI  1.2  mechanism. 

Fig.  2  shows  a  qualitative  view  of  the  pathway 
of  methane  conversion  to  C02  in  the  model.  Only  a 
very  small  fraction  of  the  CH30  is  converted  to 
methanol  by  reaction  with  water,  so  this  path  does 
not  appear  in  the  figure.  In  addition,  because 
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methanol  is  significantly  more  reactive  than 
methane,  it  is  not  accumulated.  In  fact,  no 
significant  amount  of  transient  non-radical  species 
are  accumulated,  other  than  CO  and  C02.  This  is  in 
contrast  to  the  results  for  methanol  oxidation  where 
formaldehyde  and  hydrogen  peroxide  are  produced 
in  significant  quantities. 

Fig.  3  shows  a  more  complicated  flux  diagram 
for  the  oxidation  of  methane  using  a  mechanism  [5] 
that  includes  the  formation  of  CH302  by  the  reaction 

CH3  +  02  +M  <=>  CH302  +  M  (4) 


additional  pathways  associated  with  temperatures 
below  450°C. 

02  mole  Equivalence  keff  (s- 1  )a  Fitted  keff  (s-1  )b 
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The  route  by  which  methane  is  oxidized  is  very 
different  if  this  key  reaction,  and  the  subsequent 
reactions  of  CH302  as  an  oxidizer  with  fuel  species, 
are  included.  These  reactions  are  not  included  in  the 
GRI  mechanism. 


Figure  3.  Flux  diagram  of  major  carbon  species  in 
the  oxidation  of  CH4  at  500  °C  and  25.0 
MPa  in  supercritical  water  predicted  by 
the  mechanism  described  in  Schmitt  et  al. 


The  performance  of  the  GRI  mechanism  at  these 
conditions  can  be  evaluated  on  aspects  other  than 
fuel  consumption  rates.  The  global  fit  from  Ref.  3 
suggests  that  there  is  an  oxygen  concentration 
dependence  of  order  0.66.  Table  1  shows  GRI 
mechanism  predictions  for  varying  the  equivalence 
ratios  at  constant  fuel  mole  fraction  of  0.0005  at  600 
°C.  The  calculation  reveals  an  oxygen  concentration 
exponent  of  0.734,  in  acceptable  agreement  with  the 
data.  However,  this  02  dependence  is  preserved  at 
lower  temperature  in  contrast  to  the  experimental 
results  from  Steeper  et.  al.  In  contrast,  the  same 
calculation  using  the  Schmitt/ Alkam  model  shows 
essentially  no  dependence  on  oxygen  concentration. 
The  inclusion  of  the  rapid  equilibrium  with  in  the 
Schmitt/Alkam  model  serves  to  remove  any  oxygen 
dependence  from  the  mechanism.  This  suggests  that 
the  GRI  mechanism,  when  compared  to  the  low 
temperature  data,  is  not  fully  representing  the 


Table  1 .  Effect  of  Oxygen  Concentration  on  Rate  of 
Methane  Oxidation.  Initial  methane  mole 
fraction  =5.0  xlO'4,  T=600  °C,  P=250  bar. 
aDefmed  as  keff  =  1/t,  with  reaction  time,  t, 
at  [CH4]/[CH4]0=l/e.  bPrediction  from  fit 
to  keff  =  a(02  mole  fraction)^  ,  with  a= 
7.02  and  b=0.734 

It  is  possible  that  the  GRI  mechanism  is 
appropriate  at  higher  temperature,  but  it  appears  that 
the  methylperoxyl  chemistry  must  be  included  at 
lower  temperature.  This  may  be  the  origin  of  the 
slightly  higher  rates  calculated  at  the  lowest 
temperatures  in  Fig.  1. 


Figure  4  Comparison  of  experimental  results  from 
of  methanol  oxidation  at  450  °C  and  500 
°C  at  27  MPa  with  the  predictions  from 
the  GRI  1 .2  mechanism  with  and  without 
the  addition  of  Reactions  5  and  6.  The 
unmodified  calculation  shows  essentially 
no  reaction  at  450  °C. 

Recently  our  lab  has  measured  the  oxidation  of 
methanol  and  the  simultaneous  production  of 
formaldehyde  at  440  -  500  °C  [12].  Fig.  4  shows  a 
comparison  of  the  rate  of  methanol  disappearance 
predicted  from  the  GRI  mechanism  and  this 
experimental  data.  The  agreement  between  the  GRI 
mechanism  and  the  measurements  is  poor.  This  is 
primarily  because  the  mechanism  lacks  the  key 
reactions 
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CH3OH  +  H02  =>  CH2OH  +  H202  (5) 

and 

CH3OH  +  02  =>  CH2OH  +  H02  (6) 

The  role  of  Reaction  6  is  not  particularly 
significant  at  hydrothermal  conditions,  but  without 
Reaction  5,  an  unreasonable  induction  period  is 
calculated.  Fig.  4  also  shows  the  results  from  a 
calculation  when  these  steps,  and  reverse  reactions, 
are  added  to  the  scheme  with  the  following 
parameters  used  for  Reaction  5:  A f=  3.98xl013,  Eaf 
=  19.4  kcal/mol,  bf  =  0.00,  Ar  =  3.13xl015, 
Ea  =  10.75  kcal/mole,  br=  -0.90;  where  f  and  r  refer 
to  the  forward  and  reverse  reactions,  and 
k=AexpTb(Ea/RT).  Agreement  with  the  data  is 
greatly  improved,  although  does  appear  to  be  still 
about  two  times  slower.  Note  however,  that  at  these 
conditions,  this  difference  in  rate  corresponds  to  a 
temperature  difference  of  about  10  °C.  The  results 
from  the  calculation  also  show  good  quantitative 
agreement  for  th4  formation  of  formaldehyde  as  an 
intermediate  at  concentrations  up  to  18%  of  the 
initial  methanol  mole  fraction. 

Conclusions 

The  recently  developed  GRI  1.2  methane 
oxidation  mechanism  has  been  applied  to  the 
oxidation  of  methane  and  methanol  by  oxygen  in 
water  at  25  MPa  and  temperatures  ranging  from 
400  -  630°C.  These  conditions  deviate  substantially 
from  the  much  higher  temperature  and  lower 
pressure  conditions  at  which  it  was  designed  and 
optimized.  The  results  for  the  oxidation  of  methane 
with  no  modification  of  the  mechanism  agree  very 
well  with  the  available  experimental  results  at  these 
conditions.  However,  there  is  some  evidence  that 
the  mechanism  may  be  incomplete  when  applied  to 
the  <450°C  range;  perhaps  because  the  GRI 
mechanism  does  not  contain  CH302  chemistry.  To 
represent  properly  the  oxidation  of  methanol  and  the 
formation  of  formaldehyde,  the  simple  addition  of 
one  reaction  to  the  mechanism  was  required.  With 
this  modification,  good  agreement  was  achieved  in 
the  temperature  range  of  440°C  to  500  °C. 
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